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Abstract 

We develop a simple systematic method, valid for all strengths of disorder, to 
obtain analytically the full distribution of conductances P{g) for a quasi one di- 
mensional wire within the model of non-interacting fermions. The method has been 
used in [1,2,3] to predict sharp features in P(g) near g = 1 and the existence of non- 
analyticity in the conductance distribution in the insulating and crossover regimes, 
as well as to show how P{g) changes from Gaussian to log-normal behavior as the 
disorder strength is increased. Here we provide many details of the method, includ- 
ing intermediate results that offer much insight into the nature of the solutions. 
In addition, we show within the same framework that while for metals P(g) is a 
Gaussian around (g) 3> 1, there exists a log- normal tail for g <C 1, consistent with 
earlier field theory calculations. We also obtain several other results that compare 
very well with available exact results in the metallic and insulating regimes. 



1 Introduction 

The absence of self-averaging in mesoscopic disordered conductors leads to 
important fluctuation effects in the electronic transport properties [4]. In the 
metallic regime, this gives rise to the universal conductance fluctuations [5]. 
In this regime, for a given strength of macroscopic disorder, the probability 
distribution P(g) of the dimensionless conductance g (in units of e 2 /h), for all 
possible microscopic distributions of randomness, is a Gaussian with a univer- 
sal variance i.e. its value does not depend on microscopic details of the sample 
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but depends on the symmetry of the system only. With sufficiently large dis- 
order however, the moments of the conductance fluctuations can become of 
the same order of magnitude as the average conductance, and therefore the 
average value becomes insufficient in describing the statistical properties of 
the conductor. In the extreme case of the deeply insulating regime, P(g) be- 
comes log-normal, which means that it has a very long tail so that its mean 
value is very different from the most probable value. Thus for all strengths of 
disorder beyond the metallic regime, one must consider the full distribution of 
conductances. In particular, a very broad distribution may change the qual- 
itative nature of the Anderson transition, the metal to insulator transition 
at zero temperature even in the absence of electron-electron interactions [6]. 
The numerically obtained P(g) at critical disorder for the three dimensional 
(3D) Anderson transition is indeed very broad and highly asymmetric [7]. The 
conductance at the integer quantum Hall transition is also expected to have 
a very broad, almost flat, distribution [8,9]. 

There is no theoretical method currently available to obtain directly the full 
distribution of conductances for all strengths of disorder in arbitrary dimen- 
sions [10]. Moments of the distribution have been calculated in dimension 
D = 2 + e, e < 1 [11] within the supersymmetric non-linear sigma model 
framework [12] but the critical distribution obtained from the moments [13] 
in D = 2 + e differ qualitatively from the numerically obtained P{g) in 3D. 
In the present work we will consider the distribution of conductances for the 
simpler case of a quasi one dimensional (quasi ID) wire, for which the width 
W is smaller than the elastic mean free path /, and much smaller than its 
length L. Although such a system has no phase transition, it has well defined 
metallic and insulating regimes, and a smooth crossover between them. The 
distribution in the crossover regime, where the localization length is of the 
order of the system size, should give a reasonably good qualitative picture of 
the distribution near the critical regime in higher dimensions, where the local- 
ization length diverges and becomes of the order of the system size. Even for 
this simple quasi ID geometry, only the first two moments of P(g) have been 
obtained for all strengths of disorder [14], using the non-linear sigma model. 
Numerical [9,15] as well as experimental [16] results show a broad asymmetric 
distribution qualitatively similar to 3D. 

We will use a scattering approach in which the conductance is given in terms 
of the transmission eigenvalues Tj of a finite conductor sandwiched between 
two ideal leads [17]. The finite width of the lead quantizes the perpendicular 
momenta into N values, providing N channels of scattering and therefore 
N number of transmission levels. The dimensionless conductance g is simply 
given by the total transmission probability 



N 



9 = T,Ti- 




i=i 



2 



The statistics of these levels as a function of the length L ^> I of the conductor 
is described by the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation [18]. 
This equation has been shown [19] to be equivalent to the non-linear sigma 
model [12] obtained from the microscopic tight binding Anderson model, where 
electrons hop to neighboring sites that have random energies chosen from a 
uniform distribution. An advantage of using the DMPK equation is that its 
solution i.e. the joint probability distribution (jpd) of all the transmission 
eigenvalues is known [20] and since the conductance is simply related to these 
eigenvalues, it is possible to write down the full distribution of conductances. 
The major problem is that the jpd involves iV constrained integrals for the N 
eigenvalues, and that we would be interested in the case of a large number of 
transmission eigenvalues. 

We have developed a simple systematic method, valid for all strengths of dis- 
order, to evaluate analytically the full distribution of conductance for a quasi 
Id wire, starting from the solution of the DMPK equation. For simplicity, 
we restrict our calculation mostly to the unitary symmetry case where time 
reversal symmetry is broken, e.g. by the application of a magnetic field. The 
method is valid for other symmetry classes also, and we will discuss extension 
of the calculation to these cases as well. As a check of the essential framework, 
we will show analytically that the method reproduces the known exact results 
for the mean, variance as well as the distribution of the conductances in both 
the metallic and insulating limits. We will also compute the leading correc- 
tion to the mean conductance in the metallic regime, which agrees very well 
with the exact value obtained in [14] . Since the leading order correction to the 
variance vanishes for the unitary symmetry class, we will evaluate it for the 
other two symmetry classes; again the results agree very well with the exact 
results. Also, it has already been shown in [2] that the mean and the variance 
as a function of disorder obtained within the method agrees qualitatively well 
with exact results for all strengths of disorder including the crossover regime. 

Some of the novel features of the distribution obtained within our method 
have been reported earlier [1,2,3]. Here we discuss details which provide in- 
sights about the rare but allowed configurations of the transmission eigen- 
values obtained from the method and allow us to study these features. In 
particular, we provide details that led to the prediction of a 'one-sided' log- 
normal distribution near the crossover regime in [1], the changes in the shape 
of the distribution across the crossover regime obtained in [2] as well as the 
non-analyticity in the distribution near g = 1 predicted in [3]. We also show 
for the first time within DMPK that while for metals the distribution is a 
Gaussian around (g) ^> 1, there exists a long non-Gaussian tail for g <C 1, 
consistent with earlier field theoretic calculations [11]. 

The paper is organized as follows. In section 2 we review the DMPK equation 
and its general solution. In section 3 we introduce the basic strategy of our 
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Fig. 1. Quasi-one dimensional disordered wire attached to two electron reservoirs 
fii, fi2- The DMPK equation describe the evolution of the joint distribution of the 
transmission eigenvalues Tj when system size L is increased e.g. by <5L 

new approach which is used to evaluate P(g) from the jpd of the Tj for all 
strengths of disorder in terms of certain saddle point free energies. Our method 
is based on separating out an appropriate number of discrete eigenvalues and 
treating the rest as a continuum. For simplicity, we start with separating out 
only one eigenvalue in section 3. In section 4 we obtain explicit forms for the 
density of the continuum and various free energy terms. In this section we 
also discuss some of the special features in the density and their implications 
for the distribution. In section 5 we focus on the metallic regime and obtain 
several analytical results and compare them with known exact results. We 
also obtain the non-Gaussian tail of the distribution. In section 6 we obtain 
approximate analytic expressions for the distribution in the insulating and 
crossover regimes based on the separation of one eigenvalue. In section 7 we 
extend our formulation by separating out two eigenvalues, which have been 
used to obtain the non-analyticity in P(g) near g — 1 and also to evaluate 
the distribution across the crossover regime. In section 8 we briefly discuss the 
case of other symmetries, and obtain corrections to the mean and variance of 
g and compare them with known exact results. We provide a summary and 
discuss conclusions in section 9. 



2 Brief review of DMPK equation 



For a quasi ID wire with N channels, as shown in Fig. 1, the probability 
distribution p(X) of the variables \ = (1— T^/Tj, where Tj are the transmission 
eigenvalues, satisfies the well known DMPK equation [18] 



dp_ 

dL 
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N+lJ(X)^dX a 



X a (l + X a ).J(X) 



dp(X) 



(2) 



where / is the mean free path, 
depends on the symmetry of the 
with time reversal invariance, f3 
symmetry is broken by e.g. the 
symplectic symmetry where e.g 
equation has been solved for all 
is complicated, but simplify in 



■J(X) = Ui<j \\ — Xj] 13 and the parameter (3 
ensemble [24] : (3 = 1 for orthogonal symmetry 
= 2 for unitary symmetry where time reversal 
application of a magnetic field, and /3 — 4 for 
. spin-orbit scattering is present. The DMPK 
three values of (5 [20,21]. The general solution 
the metallic and insulating regimes. In these 
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two regimes it can be written in the general form 

p ( A ) = _Lexp[-/3tf(A)], (3) 

where Z = f Ilj d\ exp[— (3H] is a normalizing factor independent of Xi, and 
H(X) may be interpreted as the Hamiltonian function of N classical charges 
at positions Xi, interacting with a two-body repulsive interaction but confined 
by a one-body confinement potential, given as 

N N 

H = J2u(X l ,X J ) + J2V(X l ). (4) 

i<j i 

We will use the variables Xi = sinh 2 Xj, in terms of which the solution in the 
metallic regime 1 <C L/l <C N is given by 



1 1 

u(xi, Xj) = — In I sinh 2 Xi — sinh 2 Xj\ In \xj — x?\, (5) 

7/ 1 

V<kX ^ = 2143^ ~ ^ ln l x i sinh2 ^l-- 
Here 7 = (3N + 2 — (3. In the insulating regime the solution is given by 



1 1 

u ins (xi, Xj) = -- In I sinh 2 x { - sinh 2 Xj\ - — In \x\ - x 2 |, (6) 

Vixj) = —r~:xf In \xi sinh 2xi I In X;. 

v ' 2Lf3 1 2(3 1 1 2(3 

The parameters N, L and I appear in the above expressions in the combi- 
nations L/l and V = Nl/L. We will consider the limit where both N and L 
approach infinity, keeping V fixed. The limit T ^> 1 will correspond to the 
metallic regime while r < 1 will correspond to the localized regime. Note 
that for (3 = 2, the solution in the insulating regime differs from the solution 
in the metallic regime only by a logarithmic term in the one-body potential. 

In terms of the variables A« or Xi, the dimensionless conductance is given by 



+ Xi i cosh x 



In the insulating regime, (g) 1 and x^ ^> 1. In this limit both lnxj and 
In |x 2 — x 2 | terms in Eq. (6) are negligible compared to the other terms. Since 
these are the only terms that differ from Eq. (5), we can assume that the 
metallic expression (5) is valid for all regimes. 
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In the metallic regime (g) ^> 1, the Aj are very close to each other so that 
a continuum description can be used with a density of A points finite be- 
tween zero and an upper cutoff given by the normalization condition. Strictly 
speaking, Eq. (3) is valid only for x 2 > 2//3T [21], but this approximation 
works well in the metallic regime. Within a global maximum entropy ansatz, 
such a continuum approximation was used [22] to show the relevance of a 
Wigner-Dyson like random matrix ensemble [25] for the transmission levels. 
The density of levels obtained from the DMPK solution within such a contin- 
uum approximation yields the correct value of the mean and the variance in 
the metallic regime [20], known from diagrammatic perturbation theory [5]. 
In the deeply insulating regime on the other hand, all Aj are exponentially 
large and separated exponentially from each other, and conductance is dom- 
inated by the smallest eigenvalue. This approximation yields the log-normal 
distribution of the conductance in the deeply insulating regime [23,20]. It is 
clear however that none of the above descriptions can be used in the crossover 
regime, where the smallest eigenvalue is neither zero, nor exponentially large. 
Our approach combines the essential features of the two descriptions and pro- 
vides a simple and systematic procedure to study the conductance distribution 
P(g) at all strengths of disorder. 



3 Generalized saddle point approximation for P(g) 



Given the joint probability distribution p(A) and the definition of g Eq. (7), 
the distribution of conductance can be written as 



1 r°° dr r°° N \ I N 1 

p ^ = zLrJ a u, dX - ex T r?TTA" 



For later use, we define a "Free energy" 



p(X). 



(8) 



(9) 



such that the distribution of g can be written as 



P(g) 'r *> r 

Z J-oo 2n Jo 



N 



n^exp [— F(A; r) 
i=i 



(10) 



This defines a statistical mechanical problem of N classical charges in one 
dimension, at an inverse temperature (3 and described by the free energy F, 
which includes a source term proportional to r. The basic idea of our method 
is to first separate out an appropriate number of the lowest eigenvalues, and 
treat the rest as a continuum by introducing a density for those eigenvalues. 
In [1] a single eigenvalue was separated out, which was sufficient to obtain the 
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qualitative nature of P(g) in the metallic regime, in the crossover region on 
the insulating side as well as deep in the insulating regime. In [2,3], the im- 
portance of separating at least two eigenvalues in order to study the crossover 
region quantitatively was emphasized. However here we will consider in detail 
the case of separating out one eigenvalue in the metallic regime since it allows 
us to obtain many of our results analytically, which provides much insight 
into the nature of the solutions, and two eigenvalues in the crossover and in- 
sulating regimes where it becomes essential. The 'density' of the continuum, 
including the source term, is obtained following Dyson [24] by minimizing the 
free energy functional, which then yields a saddle point free energy. The qual- 
itative nature of the density gives important information about the problem 
and will not only provide guidelines for appropriate approximations, but will 
also suggest novel possibilities in the distribution. We will see that at least 
for (3 — 2, the r dependence of F(A; r) is only quadratic, so the integral over 
r in (10) can be done exactly. The evaluation of P(g) is then reduced to two 
integrals in the metallic regime and three integrals in the crossover and in- 
sulating regimes, over all possible positions of the separated discrete levels, 
as well as over the beginning of the continuum. These integrals can then be 
done approximately analytically in some limiting cases to obtain qualitative 
results, or numerically to obtain quantitative results for comparison with di- 
rect numerical simulations. While we will mainly restrict our considerations 
to the case (3 = 2 only, we would point out the results for other values of (3 
whenever they can be obtained without too much difficulty 

In the following, we will discuss the steps in more detail. 



3. 1 Separation of the lowest eigenvalue 

We start by separating out only the lowest eigenvalue from the rest. This will 
be sufficient to study the metallic regime. In Sec. 7, we will show how to add 
a second discrete level, and write down the rules to extend the one-discrete 
level formulae to the two-discrete level case. This will keep the discussion of 
the basic framework simpler. 

In the variable Xi, we get 

H = H 1 + £ ufaxJ+^Vixi), (11) 

Ki<j i>l 

where 

H^^ufaxJ+Vfa). (12) 

i>l 

The idea is to treat the rest as a continuum, beginning at x 2 > x 1 . 
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3.2 Continuum approximation for the remaining eigenvalues 



Before going to a continuum description, we note that u(xi, Xj) can be written 

as 

u(xi, Xj) = u'(xi — Xj) + u'(xi + Xj), (13) 

u'(x) = — - ln(2:r sinhrr). 

The continuum description will then include the Xi = Xj terms in u'(xi + Xj) 
which are not included in the discrete part. In order to get rid of these terms 
we rewrite the sum Eq. (11) as 



H = H 1 

.(14) 



i>i z 



u'(xi - Xj) + Y u \ x i + x j) - J2u'(2xi 



The Xi = Xj terms in u'(xi — Xj) are infinite, and will not make a difference in 
the continuum limit. Using in addition the expressions for V{xj) and -u'^Xj), 
we obtain 

H = H ^ + o E x i) + E IM^) + 2« H^i sinh 2 Xi )} , (15) 

^ i,j>l i>l 



where 



Vrixi) = \x\ a = \(l - -\ (16) 



Thus for (3 = 2, the term proportional to a drops out. 

The continuum approximation can now be made, giving for (3 = 2 : 



1-0 no 

H(x±, X2, cr(x)) — V(xi) + / dxa{x)u{x,x\) + I dxa(x)Vr(x) 

JX2 JX2 

+ - ( dx I dx a{x)u{x, x')a{x') (17) 

2 Jxo J XI 



1 r b 

1X2 JX2 

where 



r 1 

V(x x ) = - x\ - — ln(xi sinh 2x x ) (18) 

and the 'density' o(x) has to be obtained in a consistent way, subject to the 
condition that it is positive everywhere. The upper limit b in principle should 
take care of the normalization condition 

f a{x)dx = N - 1. (19) 

J X2 
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However, we will be interested in the N — > oo limit, and since the contributions 
to the conductance from large eigenvalues are negligible, we will not worry 
about enforcing the normalization condition and replace the upper limit by 
oo. The Free energy can now be written in the continuum approximation as 

F(xi, x 2 ; cr(x)) — (3H(xi, x 2 ; cr(x)) H h it / — -^—dx (20) 

cosh x\ Jx 2 cosh x 

and P(g) is given by a functional integral on the generalized density a(x): 

1 r°° d,T ■ r°° r°° r 
P (9)=r ^T elT9 dx ' dx 2 D[a(x)] C xp[-F(x l ,x 2 ;a(x))}. (21) 

Zj J-oo 271 JO Jxi J 



3. 3 Integral equation for the saddle point density 



We now minimize the free energy functional by taking a functional derivative 
with respect to a(x) keeping x\ and x 2 fixed: 



5F(x 1 ,x 2 ;a(x)) 
5a(x) 



= 0. 



(22) 



This gives an integral equation for the saddle point density. We will see later 
that because the saddle point density goes to zero close to the origin, we will 
need to consider fluctuation corrections to the saddle point function. 

Following Dyson (generalized for the non- logarithmic interaction case in [20]), 
we have the saddle point equation for a(x): 



X2 



dya(y)u(x,y) + 4a 



In a(x) + Su(x, y) + — x 2 

Lj 



V tot (x) + const, (23) 



where 5u(x, y) = u(x, y) + In \x — y\. Here V tot now includes 



V t ot{x) = Vr(x) + u(x, xi) + 



IT 



cosh x 

Again the term proportional to a drops out for (3 — 2. 



(24) 



3.4 The shift approximation 



We define 

u\ = In I sinh 2 x — sinh 2 y\; u 2 = In \x 2 — y 2 \, (25) 
so that u(x,y) = —\{u\ + u 2 ). Then Eq. (23) for (3 = 2 becomes 

POO 

/ dya(y)[u l (x,y) + u 2 (x,y)] = 2V tot (x) + const; x > x 2 . (26) 
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In the limit x 2 — > 0, the density can be extended symmetrically to negative 
values which makes it possible to invert the kernel and obtain a(x). Here the 
finite lower limit poses difficulty A shift in the variables x and y to make 
the lower limit zero does not help, because the kernels are not translationally 
invariant. However in the variable A, the kernel ui(X — A') = ln|A — A'| is 
translationally invariant while u 2 (X, A') = In |arsinh\/A — arsinhv^V| remains 
invariant for a metal (A <C 1), and is negligible in the insulating regime (A ^> 1) 
compared to U\. This suggests the following shift approximation to leading 
order, to which systematic corrections can be made: 

(i) Write the integral equation in variable A 

rOO 

/ dX'[ Ul (X - A') + u 2 (X, X')]p(X') = 2V tot (X). (27) 
J A2 

where p(X')dX' = a(x)dx. Now shift to 77 = A — A2, and 7/ = X' — A2 to make 
the lower limit zero: 

roo 

/ dr)'[u 1 {r)-r)')+u 2 (r ] + X 2 ,r)' + X 2 )}p{r)' + X 2 ) = 2V tot (r) + X 2 ) . (28) 
Jo 



[ii) Write u 2 (i] + A 2 , i]' + A 2 ) = u 2 {i], i]') + Au 2 where 



Au 2 (r),r)'] X 2 ) = In 



arsinh a/?7 + X 2 — arsinh y/rf + X 2 
arsinh 2 ^^ — arsinh 2 y/rf 



(29) 



For both r\ ^> X 2 and r\ <^ X 2 , Au 2 is negligible. In the metallic regime, the 
leading term of Au 2 is a small correction linear in X 2 . So, for X 2 <^ 1 it is 
given by 



Au 2 (7],r]';X 2 ) 



A 5 



arsinh ^/fj — arsinh \fr( 



arsinh ^/fj arsinh y/rf 



?7(l + 77) Jr/(l + rf) 



■ (30) 



In the opposite limit i.e. in the insulating regime (X 2 ^> 1), the entire term 
u 2 is negligible compared to u\. Therefore we will keep the leading term of 
A«2 as a small correction to the free energy in the metallic regime only. Its 
contribution to the saddle point density itself would be negligible. In this 
spirit, we make the approximation: 

roo roo 

/ dr,'Au 2 {ri, r/; A 2 )p(t/ + A 2 ) w / d V 'Au 2 (ri, 77'; X 2 )p sp (r]' + A 2 ), (31) 
Jo Jo 

where the saddle point density p sp is the density obtained from the saddle 
point integral equation in the absence of the Au 2 term. Then the integral 
equation (27) can be rewritten with the Au 2 term taken to the right hand 
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side to add to the V to t, giving 

roo 

/ dr ] '[u 1 (r ] -r ] ')+u 2 (r ] ,r ] , )}p sp (r ] ' + X 2 ) = 2V eff ( m X 2 ), (32) 
J o 

where the effective potential V e ff is defined as V e ff(r); A 2 ) = V to t(jl + A 2 ) — 
Viiv'i ^2), w ith Hot 0*0 defined in (24) and 

V 2 (??; A 2 ) = -J Q d V 'Au 2 { v , ?/; X 2 )p sp (v' + A 2 ). (33) 

In the crossover and insulating regimes, we will take V 2 = 0, since Au 2 is 
negligible in these regimes. 

(iii) Now change the variables in the saddle point equation (32): 

sinh 2 1 = T) = X — X 2 , sinh 2 s = r{ = X' — X 2 (34) 



/ ds [in I sinh 2 1 - sinh 2 s\ + In \t 2 - s 2 \\ a(s) = 2V eff (r}(t) + A 2 ). (35) 

J 

where we have defined 

dr( 

a(s) = p sp (r]' + ^2)-^ = 2 sinh s cosh s p sp {rf + A 2 ). (36) 

Defining a(— s) = a(s) we extend the density symmetrically to negative s and 
using the relation: sinh 2 1 — sinh 2 s = sinh(t — s) sinh(t + s), the integral (35) 
can also be extended to —00, giving 

/oo 
ds[\n I sinh(t -s)\ + In \t - s\]a(s) = 2V eff (rj(t) + A 2 ). (37) 

Since r)(—s) = rj(s), this is true for all s. 

(iv) Taking a derivative on both sides of the integral equation, we get 

/oo A 
^ dsK(t - s)a(s) = 2-V eff ( V (t) + A 2 ) (38) 



where the kernel 

Kit - s) = coth(t - s) + — !— . (39) 

t — s 

This kernel K(t) can be inverted to give the saddle point density 

/oo A 
dsK-\t - S )—V eff (r,(s) + A 2 ). (40) 
-00 ds 
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Using Fourier transform, we invert the kernels K\{t) = 1/t and K 2 {t) = cotht 
to obtain Ki(q) = msign(q) and K 2 (q) = in coth(7rg/2). Then using a(q) = 
[K^q) + Kiiq)]- 1 2V^ ff (q) we write 

/oo An p—iqt 1 ^oo 1 
^-—^-—— = -/ dqsm(qt)Ul-e^ q ). (41) 
-oo 2irK 1 (q) + K 2 (q) tt 2 Jo H yH ' 2 y ' y ' 

Finally, introducing Eq. (41) into (40), the saddle point density is given by 

2 r°° r°° d 

<r(t) = — dqcos(qt)(l-e-^) dssm(qs)—V eff ( V (s) + X 2 ), (42) 

7T Z JO JO d.S 

using the fact that V e // is even in s. Note that for (3^2, there would be 
additional terms in V e ff proportional to the parameter a. 



3.5 Contributions to the Free energy 



From the saddle point density, we obtain the saddle point Free energy by using 
the right hand side of the integral equation (23) to eliminate one of the two 
integrals of the last term in Eq. (17). Substituting this result into Eq. (20) 
and changing to the A variable (since shift has to be made in the A variable), 
we write the Free energy as 

3 r°° 

F sp (Xi, A 2 ) = £ / d\V tot (X)p sp (X) + /?y(Ai). (43) 

Shifting by A2 to make the lower limit zero, and then going to the variable t, s 
and using the expression for the density, we get 



3 r°° 

F sp (X l ,X 2 ) = ^J Q dtV tot ( V (t) + X 2 )a(t)+(3V(X 1 ). (44) 



P 

fspKM, * 2 ) — 
We will write 

VtoMt) + A 2 ) =V r + V U1 + V U2 + V z (45) 
with corresponding densities a(t) = <xr + c ui + a U2 + a z (Note that for (3 7^ 2, 
Vtot has to be modified, Sec. 8), where 



V r = - (sinh 1 yj sinh 2 1 + A; 



1- / . ., 



V U1 = -- In (sinh" t + A 2 - AiJ , 
K 2 = -^hi 

V,= ,2^, ; * = (46) 
cosh t + A 2 



sinh 1 \J sinh 2 1 + A 2 J — (sinh 1 J~\~i 
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In addition, we define a density component a 2 , generated by V 2 as defined in 
(33). So the total free energy would be a sum of all combinations of partial 
free energies 

F = Y.F ab + F 1 (47) 

a,b 

where a, b runs through the four indices T, u±, u 2 , z as defined above and b runs 
through an additional index 2, and 

r°° 

F ab = ^ dtV a ( V (t) + *2)cT b (t), (48) 
2 jo 

with the partial density 

a b (t) = — / dqcos(qt)(l - e^ q ) / dssm(qs)—V b (rj(s) + A 2 ). (49) 
7T Z jo Jo as 

Here F\ = 0V(Xi). By doing a partial integration for both the t, s integrals, 
we see that the boundary terms are independent of A, giving F ab = F ba . Also, 
we note that we need to obtain the partial densities in order to calculate V 2 . 



3.6 The action for P(g) 



Since V e ff and therefore a are both linear in r for = 2, the free energy is 
quadratic in r and can be written in the form 



F sp = F° + (tr)F' + y -^-F" (50) 
and the integral over r in (21) can be done exactly. The result is 

1 /"CO POO 

P(g) = - d\J d\ 2 e- s , (51) 

Z JO JAi 

where the saddle point action S is given by 

S = -^(9-F') 2 + F°. (52) 



Note that the quadratic form in (50) is no longer true for f3 7^ 2, because of 
the term m(cx(j;)) in (23) which contains r inside the logarithm. However, we 
will see later that the leading terms in the free energy are still quadratic in r, 
and will allow us to calculate the leading corrections to (g) and var(g). 
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3. 7 Mean and variance of g as function of T 



We can obtain the moments of P(g) directly from the free energy instead of 
obtaining the full P(g) first. The nth moment of conductance is given by 

(g») = / dgg n P{g)/ / dgP(g), (53) 

JO JO 

defining the integrals as 



h(xi,x 2 ) = I dgg k e (54) 



which can be done analytically: I = y/2\F"\f (Q); h = \F"\e- Q + F'I ; 
h = ~F' 2 I + 2F'h + (2\F"\f/ 2 f 2 (Q), where f Q (Q) = (0F/2)[l + Erf(Q)] with 



Erf(Q) the Error function, f 2 (Q) = (f (Q) - Qe~ Q )/2, and Q = F'/y/2\F"\ 
. Defining 

J k = dxi / dx 2 e~ F h(xi,x 2 ), (55) 

JO Jx 2 

we have that the average and variance of g can be expressed as 

(g) = Jx/ J ; var(g) = J 2 /J - (Ji/J ) 2 . (56) 



4 Density and Free energy 



In order to evaluate the free energy terms explicitly, we first need the five 
components of the total density, given by Eq. (49). While a z and a Ul can be 
evaluated exactly analytically, the others require approximations. In particu- 
lar o 2 can be evaluated only after the potential V 2 has been obtained from the 
other partial densities, and V 2 itself involves approximations. In principle, one 
could evaluate all the free energy terms directly numerically without any fur- 
ther approximations. This will involve at most triple integrals, Eqs. (48) and 
(49). More importantly, we will see that the qualitative nature of the various 
partial densities provide much insight into the nature of the problem and not 
only provides guidelines for appropriate approximations in various regimes, 
but also suggests novel possibilities in the distribution. We therefore will eval- 
uate the density integrals analytically, albeit approximately when necessary. 
Here we list the results for the different cr's; details of the calculations are 
given in Appendix A. 

For a z (t) and cr Ul (t) the results are exact while for cr U2 (t) and cr(t) good 
approximations can be obtained analytically by modelling the derivatives of 
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the potentials Vr and V U2 as 

V£(s) « Ts 



04 



cosh 2s + cosh Si 



(57) 



where a\ and 5i are chosen to match the exact form at s — > and at s = x 2 , 
and 



s 2 + si 



as: 



s 2 + si 



(58) 



(s 2 + s^) 2 + 4cV 

with a and c obtained by matching (s) at s = and s = s m , s m being 
chosen to follow the maxima which occurs at s = Xd for Xd — > and at s ^ 12 
for a; 2 , >> 1, where x\ = x\ — x\. Using these approximations we list the partial 
densities: 



<Tz(t) = 



x 2 + t 



+ 



x 2 -t 



v/a 2 (i + a 2 ) [((x 2 + ty + (n/2yy ((x 2 - ty + (tt/2) 2 ) 2 



<? U1 (t) 



1 

7T 



(tt - oo)/2 



+ 



(7r + a )/2 



(tt - a ) 2 /4 + t 2 (7r + a ) 2 /4 + t 2 



(59a) 
(59b) 



<?u 2 (t) 



Pi + c 



TT + Pi + C 



2tt/3i + (A + C) 2 



t 2 + (7T + A+ C ) 2 / 
7T + ft - C 



27T/3i \t 2 + (Pi - c) 2 * 2 + (TT + Pi - cf 



(59c) 



a T {t) ~ r 1 + 



04(^1 



7T 2 + (5 2 



04 



4 sinh Si (t 2 + (tt 2 - 5i)/4) 2 + (tt^/2) 2 cosh 2t + cosh Si J 

(59d) 

with the following parameters: 



a = cos _1 [2(A 2 - Ai) - 1]; A± = Pi ± c + as 2 m /2c; a=l 



X2P2 



2 2 2 i 1 Slllll / o 2\ 

Pi = c + s m = 2:2 smh p 2 ; p 2 = ^ {x 2 - ay. 



^.i Positivity of the density 



Both <j ui and <t U2 diverge with negative coefficients in the limit t — (i.e. 
A = A 2 ) and A 2 — Ai — > 0. This means that the total density (without the 
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t 1 1 1 1 1 1 1 1 1 1 n 4^ i i i i i i i i 

0.2 0.4 0.6 0.8 1 1 2 3 4 5 6 7 

t t 

Fig. 2. Total density atot(t) as a function of t (a) in the metallic regime T = 5 
with fixed x\ = 0.1 for three different values of x 2 and (b) in the insulating regime 
r = 1/3 with fixed x\ = 0.5, for three different values of X2, showing that the 
negative part of the density gradually disappears with increasing x 2 . Note that in 
the metallic regime, a(t) saturates at T for large t. 

source term a z ) given by 

cr to t(t) = o- r + a ui + a U2 (60) 

will always become negative near t — in this limit. Negative total density 
for any t is clearly an unacceptable solution. To determine what to do with 
these negative density solutions, we need to understand why our saddle point 
equations produce such solutions. Both a Ul and a U2 are contributions from the 
repulsive interaction terms. The total density becomes negative when A 2 — 
Ai becomes less than certain value depending on T. Thus the negative total 
density corresponds to the unacceptable configurations of A 2 too close to Ai, 
which are very costly in energy due to the mutual repulsion of the eigenvalues. 
It is therefore clear that such solutions should be discarded. In fact, positivity 
of the total density can be guaranteed simply by imposing a minimum on the 
allowed values of A 2 — Ai, such that a tat {t = 0) > 0. In the metallic limit 
A 2 <C 1, and the condition of positivity becomes 

x 2 2 -x\>(^)\ (61) 

while in the insulating limit A 2 > Ai > 1 we have 

x 2 -x 1 >l/2T. (62) 

Figure 2a shows the total density cr tot (t) for x 1 = 0.1 and three different values 
of x 2 in the metallic regime; it shows how the negative part of the density for 
small t gradually disappears with increasing x 2 , such that the density can 
always be made positive for all t by choosing the appropriate lower cut off x r 
for x 2 . In all our analytical evaluations on the metallic side we will use the lower 
cut off (61). Note that since x 2 — X\ 1 in the metallic limit, the restriction 
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Fig. 3. atot(t = 0) as a function of %2 for two different values of T: (a) in the metallic 
regime, for a given fixed x± = 0.1 and (b) in the insulating regime, for a given fixed 
X\ = 1.0, and two different values of T, showing how the cutoff for the positivity of 
density depends on disorder. 

(61) also imposes a condition T ^> 2/tt for which a positive density metallic 
solution can be obtained. This gives a crude estimate that the metallic region 
ends and the crossover region begins at T as 1/2. Similarly, Fig. 2b shows 
<Jtot(t) for x\ = 0.5 and three different values of X2 in the insulating regime. 
In figure 3a, a(t = 0) is plotted as a function of X2 for two different values of 
T in the metallic regime for a fixed x\ = 0.1, while Fig. 3b shows the same 
in the insulating regime for a fixed X\ = 1.0. From this figure 3b we can see 
how the increase in disorder (decrease in F) increases the cutoff for x 2 that 
ensures the positivity of the density. Fig. 4 shows several numerically evaluated 
roots X2 — x r of a(t — 0). For comparison with Eq. (61), we plot the scaled 
quantities (Tx r ) 2 as a function of (Txi) 2 which, according to Eq. (61), should 
be a straight line in the metallic regime. Figure 4 shows that the metallic limit 
agrees very well with Eq. (61), and that the same equation provides a good 
estimate and an upper bound in the crossover and insulating regimes as well. 
Nevertheless, in our numerical evaluations, for each x±, we will find the value 
of x r numerically by evaluating the total density at t = 0. Because of the 
cutoff in x 2 at x r , P(g) will be given by 



Note that since we used the variable x instead of A, the action S(xi,X2) will 
contain appropriate Jacobian terms. 

The inequalities (61, 62) show that the separation between x\ and x 2 increases 
with decreasing T. Since T ^> 1 describes metals and T 1 describes insula- 
tors, this gives rise to the well-known qualitative picture that all eigenvalues 
in metals are close to the origin while they are all large and well-separated 
in the insulating regime. In addition, Eqs. (61) and (62) suggest two more 
qualitatively distinct regions involving rare but allowed configurations: 




(63) 
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)l , 1 , 1 , 1 , 1 

12 3 4 

( r*/ 



Fig. 4. Values of x 2 = x r for which atot(t = 0) scaled by T plotted as a function 
of (Txi) 2 for several T in the metallic, insulating and crossover regimes. Eq. (61) 
fits the metallic regime well, and is a good estimate and an upper bound for the 
insulating and crossover regions. 

(i) T >> 1, x± >> 1, x 2 1, but x 2 ~ Xi. We will show in this paper that even 
for metals, this gives rise to the existence of long tails in P(g) in the region 
g < 1 (Sec. 5). 

(ii) r <C 1, xi <C 1, x 2 1. In this insulating region a Gaussian cut off in the 
log- normal conductance distribution for g > 1 was found in [1] (Sec. 6). 

Thus the qualitative features of the density provide us with much insight, 
not only about various approximations and their limitations, but also sug- 
gesting rare but allowed configurations giving rise to novel possibilities in the 
distribution. 



4-2 Fluctuation correction to the functional integral 



As we have seen, the generalized density a(t) goes to zero at t — 0. This means 
that the saddle point evaluation of the functional integration over the density 
requires corrections. The fluctuation correction will involve the integral 



exp 



dt / ds5a(t)u(t, s)8a(s) 
Jo 



(64) 



where u(s,t) includes the two body terms, and 5a 
general of the form I a 



a — a sp . The result is in 



l/ydet\u(t,s)\. This will give rise to a contribution 
to the free energy of the form 5F = +| ln(det\u\). Although it should be pos- 
sible to evaluate the determinant numerically, we can already argue that the 
correction is usually negligible (see Appendix C). The fluctuation correction 
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to the free energy is relevant only in the extreme metallic regime and it is 
given by 

5F^Un\ d . (65) 



4-3 The partial free energies 



The partial free energies related with a z and a u \ given by Eq. (48) are calcu- 
lated exactly for any regime. In particular, F zz is relevant for the calculation 
of the exact variance in the metallic regime as well as to see the qualitative 
change in the solutions for large and small x 2 \ 



F = 

■*■ 7.7. 



Z 2 



1 1 



(3 sinh^(2x 2 ) L 3 Ax l + sinh 2 2rr 2 
There are also two direct terms 



(66) 



*i = f r *i ~ \ M*i smh(2 Xl )); F; = — ^— , (67) 
2 2 cosh X\ 

and the contribution from the fluctuation correction to the functional integral 
SF (65) valid for x 2 <C 1. Other free energy terms can be calculated within the 
same approximations as used for the densities (see Appendix B for details). In 
particular, the free energies can be obtained analytically in the metallic and 
insulating regimes, keeping the leading orders in x 2 . On the other hand, we 
will see later that an interesting crossover regime on the insulating side will 
be determined by the case where x 2 ^> 1, but x\ can be arbitrary. Free energy 
expressions in this regime can also be obtained approximately analytically. 
Note also that these expressions will be valid for (3 ^ 2 as well, for which 
there will be additional Vj terms and therefore additional free energy terms. 
For quantitative results on P(g), we will integrate numerically Eq. (48) to 
obtain the free energy. 

In the final saddle point action, we have written the total free energy in the 
form given in (50). Remembering a factor of two for symmetrical terms F^ = 
F ba for b 7^ a, we then have 



F — Fyt + + + F° + 5F + 2[F uir + ], (68) 

F' = -Ff + - [F Tz + F ZU1 + F ZU2 + F z2 + Fi T ] , (69) 

z z 

F" = - 2 F ZZ . (70) 

z A 

Given the free energies, the action S determines the full distribution P(g). We 

will discuss several quantitative results in the metallic and deeply insulating 
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regimes as well as qualitative results in the insulating side of the crossover 
regime analytically. The distribution in the crossover regime has to be obtained 
numerically. 



5 The metallic regime 



The metallic regime is characterized by x 2 <C 1 and the Free energy terms 
in this regime can be obtained analytically by expanding the potentials (46). 
Details of the calculations are given In Appendix B. 



5.1 Free energy in the metallic limit: 



Using x\ < x 2 <C 1 and therefore a/A^z ~ \jx\ — x 2 , the total free energy in 
the metallic limit is given by 



F° ~ aiY 2 x\ + a 2 \n(x 2 2 - x\) - a^Y^Jx 2 2 -x\- lnxi, (71) 
F'^T-hTxl + b^xl-x 2 - F"n-±;, (72) 

where a x = 3tt 2 /8, a 2 = 3/2, a 3 = 2vr, b x = 2/3 + 2 (8 - 2((2) - 3((3)) /tt 2 
and 6 2 = 32/7T 3 , where ((n) is the Riemann zeta function. 



5.2 Saddle point analysis in the deep metallic limit 



Since 12 C 1, and T ^> 1, we will keep only the dominant first term in F' 
for the purpose of obtaining a saddle point solution (we will see that x 2 is 
0(l/r), so that the neglected terms are smaller by a factor l/Y 2 compared 

to the leading term ). In terms of the scaled variables /x = Y\J x 2 — x\ and 
v = Tx\, the free energy in this limit is 

F° ss ai(/i 2 + v 2 ) + 2a 2 In jj, - a 3 /i - In v - ln/j + - ln(/i 2 + z/ 2 ), (73) 



F' w r; F" = - — . (74) 
15 V ; 

The last two terms in F° comes from the Jacobian of transformation dxo = 



[id[ij '(rV/j 2 + v' 1 ). Using (52), the solutions to the saddle point equations 
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dS/ dji = dS/ dv — are given by 



}_ = ^1±1, 1 = _ 2ai + J_. ( 75 ) 

Eliminating v we obtain a cubic equation for /x, whose real solution turns out to 
be less than the cutoff due to the positivity of the density. This means that the 
metallic solution is dominated by the boundary value \i = T^Jx 2 — x\ = 2/ir. 

The corresponding v = fi/y/2. Thus the saddle point analysis gives both x\ 
and x 2 to be of order 1/T. 



5.3 Correction to (g) in the metallic limit 



The nth moment of P(g) can be obtained from Eqs. (54) and (55). The average 
conductance in the metallic limit is T with leading correction ~ 1 /T, for /3 — 2, 
with a coefficient obtained exactly in [14]. In order to compare with [14], we 
write (g) as 



(g) 



dx x j£ dx 2 / °° dg[(g - F') + F']e 



(76) 



o dxi J Xr dx 2 Jo dge b 
The first integral over g in the numerator can be done simply and we write 



roc 

/ dg[(g - F') + F'] 
Jo 



\F"\e 



-F° 



+ F 



roo 

Jo 



(g--F') 2 pO 



2\F"\ 



(77) 



In the metallic regime, F' ~ T > 1, and = 1/15 <S 1; the first term in (77) 
is exponentially small and we neglect it. The numerator and the denominator 
in (76) then have the same integral over g. For F ~ T » 1, this integral 
depends on |F"| only. Since |F"| is a constant, it cancels out from the ratio 



and we have 



(g) 



Io° dx i $™ dx 2 F'e 



-F° 



(78) 



JO UX 1 Jx r UX 2 t. 

The leading term in F' is T, independent of x±, x 2 , so the leading term for the 
average conductance is simply T. Scaling Xi = XiT, we write 



<<?> = r 



(79) 



where 



Jo°° dx i !s T dx 2 



7/ 



b\x 2 2 - b 2 \l x^ - x\ 



-F° 



/o°° dx i !T r dx 2e 



-F° 



F° = a\x\ + a 2 ln(xl — x\) — a 3 \f x 2 



x 1 



lnxi, 



(80) 
(81) 
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and x r = yjx\ + (2/tc) 2 . Note that while the free energy expressions are valid 
only for x ; < 1, the scaled variables be very large. However, the results 

become insensitive to the upper limit beyond 4. Setting the upper limit at 4 
and doing the integral numerically we get 

r] = .027 (82) 

compared to the exact result r] = 1/45 ~ .022. This shows that our results 
are good in the metallic limit up to order 1/r. This is consistent with the fact 
that Xi ~ 1/r, (i = 1, 2), and that we kept only terms up to order X{ and Txf 
in Eqs. (71, 72). 

5.4 Variance in the metallic limit 

Using g 2 = (g-F') 2 +2gF'-F' 2 , and the above approximation that F' 2 /2\F"\ > 
1 with F" constant, we can write the second moment, following the same pro- 
cedure as above, as (g 2 ) = {\F"\ + F' 2 ). Then the variance is simply 

var(g) = (g 2 ) - ((g)) 2 = \F"\ = -L (83) 

which is the correct variance in the metallic limit. We do not calculate the 
correction to the variance for (3 = 2, because the leading correction is of order 
1/r 2 , and our approximate free energy expressions are not good up to this 
order in the metallic regime. However, the leading correction for (3 7^ 2 is of 
order 1/r, and we will see later that our approximations yield values very 
close to the exact result. 

More generally, we can evaluate the integrals Jk defined in (55) numerically 
and obtain the average and variance of g as a function of the disorder parame- 
ter T. Results of these calculations have been reported in [2]. They are in very 
good agreement with the exact results in the metallic region. 

5.5 Conductance distribution in the metallic limit 

The saddle point analysis of the metallic regime gives both x\ and x 2 of order 
1/r, which justifies keeping only the leading term in F' in Eq. (72), namely 
F' = r, deep in the metallic regime. Then the distribution P(g) becomes 
trivial, because the (^-dependent part can be taken out of the integrals over 
x\ and x 2 and the integrals only give the normalization factor. We obtain the 
known Gaussian distribution: 

P(g) oce-^-r) 2 ; (84) 
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Saddle point 



T = 2 




Fig. 5. Gaussian distribution of the conductance in the metallic regime numerically 
obtained from Eqs. (71) and (72) for T = 2, compared with the saddle point result 
Eq. (84). 

with mean (g) = T and var(g) = 1/15. Using Eqs. (71, 72), P(g) can be 
directly obtained by numerical integration. As an example, Fig. 5 shows the 
result of the numerical integration for T = 2, compared with Eq. (84) appro- 
priately normalized. 



5.6 P(g) for j<Cl m the metallic limit 



An important question for the metallic distribution is whether there are long 
(non-Gaussian) tails for g <^ 1. These will signal rare events (localized states) 
dominating the tail of the distribution, as suggested by field-theory models 
[11]. With appropriate approximations of our free energy, we can analytically 
consider this region qualitatively. We will be interested only in g <^ 1, which 
can happen only if Xi ^> 1 and x 2 1. However, metal means large T, and 
the positivity of the density requires that x\ — x\ > 4/tt 2 T 2 , Eq. (61). This 
means that x 2 can be very close to x\ 
r < 1, and x 2 > xi). 



(in contrast, for the insulating case, 



The free energy in the limit x 2 ^> 1 can be written as 

F° ~ 2r 2 ^ - 6Tx 2 + Tx\ - ]- ln(xi sinh(2xi)); 



(85) 



F' 



+ 



2T 



cosh Xi sinh2x2 



F" w — e 
3 



-4X2 



(86) 
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Fig. 6. P{g) for V = 1.0 obtained in [2] compared with normalized (88) (solid line) 
and normalized (84) (dashed line) . The log-normal tail of the distribution for g <C 1 
is clearly distinguishable from the Gaussian tail. Note that T = 1.0 is slightly beyond 
the metallic regime (a larger value would be metallic) where (84) is not strictly valid. 
Nevertheless, it clearly shows the Gaussian tail in the g 3> 1 limit as expected for a 
metal but a log-normal tail in the g <C 1 limit. 

The saddle point solution is given by 



dS _ 2(g-F') 
dx x 2 IF" I 



dF'\ dF° 

+ — = 0. 



dxi I dx\ 



*7) 



Since the quantity 1/\F"\ is exponentially large, the saddle point solution is 
given simply by putting its coefficient equal to zero, i.e. g = F' + 0(e~ Ax2 ). 



For both x\ 1 and x 2 ^> 1 but x 2 ~ X\ and T ^> 1, the second term in F' 
can become the dominant term (the first term dominates for x 2 ^> X\ ^> 1 
and r <C 1 which will describe the insulating region, and also for x\ <C 1 and 
1 which will describe the crossover region) . Then the saddle point solution 
for x 2 gives g = 2T / sinh2x 2 . The leading order solution is x 2 w |ln(4r/g), 
leading to 



P{g) oc e 



-2<9> 2 (ln 



4(g) 
9 



' 2(g) , 



((?)>!, 



where we have used the fact that in the metallic limit (g) = T. The fluctuation 
correction does not add any further g dependence. Figure 6 shows P(g) given 
by Eq. (88) compared with the distribution of g calculated from the two sep- 
arated eigenvalue case for T = 1.0, as well as the Gaussian tail expected from 
Eq. (84). Note that T = 1.0 is slightly beyond the metallic regime (a larger 
value of T would be metallic). Nevertheless it clearly shows the existence of a 
log-normal tail in the c/ <§C 1 regime compared to a Gaussian tail in the g ^> 1 
regime for the same value of T close to the metallic limit. 



24 



6 The insulating and crossover regimes within the one separated 
eigenvalue approximation 



As emphasized in [2,3], we need to isolate at least two eigenvalues in order to 
study the insulating and the crossover regimes. This we will do later. On the 
other hand, the one-eigenvalue framework allows us to obtain many analytical 
results, which helps us to better understand the qualitative nature of the 
solutions. We will therefore expect our results in this section to be valid only 
qualitatively. In particular, we will point out why the one eigenvalue separation 
fails very close to g — 1 even qualitatively. In all other cases the numerical 
evaluations with two separated eigenvalues improve the one eigenvalue results 
quantitatively, without changing them qualitatively. 

6. 1 Free energy and saddle point solution in the insulating limit 

In the limit T 1 and x 2 >> x\ >> 1, the free energy simplifies to 

F° « 2T 2 x 2 2 - 6Tx 2 + Tx\ - x 1 (89) 



F' « 4e _2xi ; F" « --e^ x \ (90) 

3 

As before, since x 2 1, the saddle point solution is given by putting g = 
F' + 0(e- 4x2 ) = Ae~ 2x \ so that x lsp = -±lnw, u = g/A. The saddle point 
value of X2 is obtained from 

|^ = 3e^(g - F'f + ^ = ^ + 0(e~^) = 0. (91) 
ox 2 ox 2 ox 2 

This gives x 2sp = 3/2T, independent of g. As we will see (Eq. 95), since 
xup = — lnu/2, Xi S p ~ l/2r, so x 2sp > x lsp for T < 1. 



6.2 Conductance distribution in the insulating limit 



Since x 2sp is independent of g, it only contributes to the normalization con- 
stant. The distribution of conductance in saddle point approximation is de- 
termined by X\ only, 

P sp (g) oc e -(-*i+ r *?) oc e -^ lnu+1 /r) 2 _ (92) 
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To this we need to add the fluctuation correction around the saddle point 



P(g) oc e~ Ssp /^\S"\, where 

2 

dx\ \F"\ \ dx\ J ' \F"\ \dxi J ' dx\ 



_0 2 S _ g-F' ( d 2 F'\ 1 fdF'Y d 2 F° 

~ f>™2 ~~ I 17// 1 I fl™2 I + I ev/ I I o„ ) o~.2 • v / 



The first term is zero at the saddle point. The third term is negligible com- 
pared to the exponentially large second term, which gives S"(x sp ) oc g 2 . The 
distribution of In g can then be written as 

P(\ng) = P{g)^- g = 9P(9) oc e^^m/r) 2 ■ (94) 

This is a log-normal distribution with mean and variance given by 

<-ln(s/4)) = i; var(\n(g/A)) = 1 (95) 

which agrees with known results. 



6.3 Mean and variance of g as function ofT 



In the insulating limit, the density is zero, and the calculation of the mean 
and the variance becomes simple. In the crossover regime, one can use the full 
free energy to numerically obtain the mean and the variance using (54-56). 
These two moments of P(g) have been reported earlier in [2] and compared 
with exact result. They fail in the crossover regime, but the two separated 
eigenvalue approximation substantially improves the agreement with exact 
results both in the insulating and in the crossover regimes. 



6.4 Distribution in the crossover regime 



To obtain quantitative results in the crossover regime, we will need to obtain 
the free energy terms by numerically evaluating Eq. (51). However, the regime 
X2 1, X\ <^ 1, corresponds to T<,1 which characterizes the crossover regime 
from the insulating side. For this case, the free energy contributions are the 
same as in Eqs. (89,90). This time, however, the fact that x 2 3> x 1 makes 
the first term in F' dominate. The results have been discussed in [1], where a 
'one-sided log-normal distribution' with sharp features in P(g) at g = 1 was 
predicted. Here as an example, we show in Fig. 7 the distribution P(g) for 
T = 5/12 calculated by using the results in [1]. On the other hand, note that 
the region T ■>!, which corresponds to the crossover regime from the metallic 
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Fig. 7. The 'one-sided log-normal distribution' of the conductance in the crossover 
regime from the insulating side, where sharp features in the distribution were found 
around g=l [1]. 

side, is characterized by x 2 ~ 1, X\ <C 1. Our approximate free energies are 
not valid in this regime. 

The saddle point solution of x\ is again given by g = F' and we have cosh x\ = 
-^=. Note that since cosh a; i > 1, the boundary of the saddle point solution is 
given by g < 1. It is the existence of this boundary that makes the current 
solution fail very close to g = 1. The two separated eigenvalue framework 
reveals the existence of a non-analyticity in the distribution near this point 
[3]. 



7 Separating out two eigenvalues 



As mentioned earlier, the two separated eigenvalue case has been treated in 
[2,3]. Here we give a summary for completeness. Separating out 2 eigenvalues 
simply means that equation (12) for Hi is now replaced by 

#1,2 = U ( X i' X l) + U ( X i> X l) + U ( X 1> X l) + V ( X l) + VM- (96) 

The continuum approximation (17) now becomes 



rb 

H(x 1 ,x 2 ,x 3 ; <j{x)) = V(xi) + V(x 2 ) + u(xi,x 2 ) + / dxa(x)u(x, xi) 

JX3 

rb | rb rb rb 

+ / dxa(x)u(x, x 2 ) -\ — dx dx'a(x)u(x, x')a{x') + / dxa(x)Vr(x), 

Jxz 2 Jx3 Jx:j JX2 

(97) 



27 



while Eq. (20) for the free energy is now 

%1~ XT /*oo (j(x) 

F(xi, x 2 , x 3 ; a(x)) = fiH(x 1 ,x 2 , x 3 ; a(x))-\ 1 Yir / — - 2 —, 

cosh x\ cosh x 2 Jx 3 cosh x 

(98) 

and the distribution P(g) becomes 

P( g ) = L / ^-e iT9 / d Xl / dx 2 / dx 3 / D[a{x)}e- F ^ x ^' u{x)) . 

Z J-oo Z7T JO Jxi Jx 2 J 

(99) 

The shift parameter in (28) is now A3 and V to t Eq. (45) has two additional 
terms 

VtoMt) + A 3 ) =V r + V un + V U21 + +V U12 + V U22 + V z (100) 

with corresponding densities a(t) = cr r + a uil + a U21 + a Ul2 + a U22 + a z . Here 
V Ull and V Ul2 are just the old V Ul and V U2 with shift parameter A 2 replaced by 
A3, (similarly for a Ull and cr Ul2 ). The new terms V U21 and V U22 can be written 
down following the substitution 

V U21 = V U1 (X 2 -> A 3 , Ax -> A 2 ); K 22 = K 2 (A 2 -> A 3 , Ai -> A 2 ), (101) 

and similarly for cr U21 and <r U22 . In terms of the new saddle point densities and 
free energies, the expression for P(g) is similar to (51), with one extra integral 
over A 3 : 

1 poo poo POO 

P(g) = - d\J d\ 2 / d\ 3 e~ s (102) 

Zi JO JAi J\ 2 

where the action S has the same form as (52). The mean and the variance can 
be obtained from the integrals Jk of (56), now given by 

roo roo roo 

J k — J dxi j dx 2 j dx 2 e~ F I k (x 1 ,x 2 ,x 3 ). (103) 

JO ill 



The extra integral over A3 makes analytical calculations of P{g) difficult. Even 
numerical evaluations become involved. We will not only put V 2 = as dis- 
cussed before, but also keep only the interaction between nearest A's, neglect- 
ing in the potential terms the argument Ai compared to A3. Thus e.g. we set 

K n ~-^m(smh 2 * + A 3 ) (104) 

Since A3 ^> Ai beyond the metallic regime, this is a valid approximation in 
the regimes we are interested in. 

The rest of the calculations are exactly as discussed for the one eigenvalue 
case. For example, the density again becomes negative for sufficiently small 
x 3 — %2- The positivity of the density is guaranteed by numerically obtaining 
the minimum value of x 3 for which the density is zero at its minimum. It is 
straightforward to repeat the evaluations for the mean and the variance of g as 
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Fig. 8. Distribution of the conductance in the crossover regime. A broad distribution 
is found for T = 1/2 (squares); eventually P(\n.g) goes to the log- normal distribution 
in the insulating regime. For T = 1 (diamonds) P(g) is closer to the Gaussian 
distribution (inset, from [2]) which is expected in the metallic regime. 



functions of T within the two separated eigenvalue framework in the crossover 
regime, and in the insulating regime where the density is put equal to zero. 
These results have been reported before [2] and shows that they are much 
better than the one eigenvalue case. In particular the variance now follows the 
non trivial hump at T = 1/2 in the exact result [14]. Details of the results 
for P(g) in the crossover regime have also been reported in [2]. They compare 
well with numerical results available [9], including the direct Monte Carlo 
evaluation of Eq. (8) [27] . As an example of the results obtained by separating 
out two eigenvalues, we show in figure 8 the distribution of the conductance 
as a function of In g and g (inset) in the crossover region. For T = 1/2, a broad 
distribution is observed, as reported in numerical studies. On the other hand, 
P(lng) goes, eventually in the insulating regime, to a log-normal distribution 
(as the one shown in Fig. 7). For T = 1, which is closer to the metallic regime, 
the profile of P(g) is close to a Gaussian distribution, as expected. 

Note that in our derivation for the distribution in the crossover region within 
the one separated eigenvalue framework, the saddle point solution existed only 
on one side of g — 1 while the other side was dominated by the boundary value. 
Separating out two eigenvalues allows us to treat the case near g — 1 with more 
accuracy. As discussed in [3], we obtain a discontinuity in P'(g) = dP(g)/dg 
at g — 1 + a where a — 1/ cosh 2 (3/2r) with P' ~ e~ 2 l v for g > 1 + a and 
with an almost flat distribution for g < 1 + a. In figure 9 we show P(g) and 
its derivative (inset) around g — 1 for two different values of disorder, both in 
the insulating regime, illustrating the discontinuity close to g — 1 reported in 
[3] and observed in numerical analysis [7]. 
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Fig. 9. P{g) and its derivative P'(g) (inset) around g = 1 in the insulating regime. 
The discontinuity in P(g) is at g = 1 + a [3]. Here a ~ 1CT 6 (1(T 8 ) for V = 1/5 

(r = i/8). 

8 The case (5 ^ 2 

In this section we come back to our simpler one separated eigenvalue approxi- 
mation. As in section 2, we will use the metallic limit Eq. (5) to be true for all 
regimes. Then the shift approximation and the corresponding expression for 
density carries through, the only difference is that there are three additional 
terms to be added to the V tot , defined in (24). One is obtained from the last 
term in (15), while considering the continuum approximation: rewriting this 
last term in the A variable, using A = sinh 2 x, and then doing the change of 
variable sinht = A — A 2 , as in Eq. (34), we have 

V«(t) = ^[^(arsinh^sinh 2 1 + A 2 ) + ln(sinh 2 1 + A 2 ) + ln(cosh 2 t + X 2 )]. (105) 

The same last term in (15) also gives a direct term from the separation of the 
lowest level, 

V^{x 1 ) = 2aln(xisinh(2xi)). (106) 

A third term comes from the 'entropy' term of Dyson [24] derived as the second 
term in (23) 

V 2 a (t) = -4a In cr(0. (107) 

It is normally assumed that this term does not change the density appreciably, 
so the density appearing in (44) can be taken to be the density due to the 
other terms. However, it includes a z , and therefore the free energy is no longer 
quadratic in z. This means that instead of doing the r integral exactly analyt- 
ically, we have to do it either numerically, or try for a saddle point solution. 
On the other hand, since the largest density in the metallic regime is a"r, we 
can expand the logarithm and keep only the leading term in z. As we will see, 
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this will be sufficient to give the leading correction to the variance of g in this 
regime. 

8.1 Correction to the metallic mean conductance 

The leading contribution to determines the leading correction to the mean 
conductance (g). We have 

Ki = f / dtVrf, (108) 

where the density is obtained from the potential V"(t), Eq. (105). Comparing 
with the potentials we already have, we can rewrite it as 

V«{t) = -2a (K 2 L^o + KJa^a, + KJa^i+a.) • (109) 

Note that the last term is a > 1 term but will contribute in the metallic 
limit when A 2 < 1. Thus the free energy can be written down in terms of 
known results 

Fzi — ~ 2« (F ZU2 \ Xl ^ + F ZUl \x d ^\ 2 + F ZUl \x d ^i+\ 2 ) ■ (HO) 
Collecting all the (already known) results, we obtain in the metallic limit 

Ki ~ az (I - ta) . (in) 

Since there is a constant term independent of Xi and this contributes to F', 
we immediately get a correction to the average conductance T, given by 

to > = r + 4 = r + (i-|)i, (ii2, 

where we have included a factor of two for the symmetric term in the free 
energy. The correction agrees with exact result. 

8.2 Correction to the metallic variance of conductance 

The leading correction to the variance comes from the contribution to F^ 
from V 2 a (t), Eq. (107). Note that a r is the largest term, except very close to 
t = in a small range where o~ ul and o~ U2 become comparable. We neglect this 
complication, and in order to extract the z 2 term, take as a first approximation 
a(t) ~ cr r (t) + a z (t). Then using o~r(t) cr z (t), we expand the logarithm, 
V^(t) ks lna r + o- z /a r . The first term to leading order will give a constant 
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which we ignore, and a term proportional to x\ in the metallic limit. The 
second term to leading order would consist of taking the A2 — > limit of both 
the densities. Then we have 

Br 1 Rn>z 2 r°° 

Ki = §/ dt(-4a) f o,; = A 

The result of this last integral is 12/n 3 . Then we have 

F « = -rw f 1 - 1) £ < 114 > 

which gives a correction to the variance 

»>=^H)^- (n5) 

The exact result [14] can be written as 5(var(g)) = -^(1 — with /i = 
8/315 = 0.0254, while for us, fi = 24/vr 6 = 0.0250, which is very close. The f3 
and T dependence are also correct (in [14], (/L = 2T; (g) and var(g) have a 
factor 2 and 4 from spin, respectively). 



7T 2 /4 - 3t 2 

(t 2 + 7T 2 /4) 3 



(113) 



8.3 The density cr" 



Since the potential Vf can be written in terms of potential already consid- 
ered, the corresponding density erf can also be written down from the known 
densities, given as 



a i(t) — — 2a (<7 U2 \ Xl ^Q + <j Ui \\ d ^\ 2 + o«i U d ^i+A 2 ) 



(116) 



In the metallic limit this simplifies to 



af(t) « 8a 



1 + 1 x 2 



4t 2 + 7r 2 2nt 2 + x 2 2 



(117) 



while in the limit x 2 — > 0, 



5 a 



+ 7<*(* - o + ) 



4t 2 + 7T 2 4 



'118) 



which agrees with 
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9 Summary and conclusion 



Given the distribution of the transmission eigenvalues, we have developed a 
simple method to obtain the full distribution of conductances P(g) at all 
strengths of disorder. The method is based on the mapping of the N eigen- 
values to an electrostatic problem of N charges on the positive line with one- 
and two-body interactions. Instead of treating all the charges as part of a con- 
tinuum in the N — > oo limit, we introduce the idea of separating out a finite 
number of charges close to the origin (corresponding to the most conducting 
channels) and treating the rest as a continuum. This allows us to treat the 
strong disorder region as well where the continuum model breaks down due to 
the large separation between the charges, and where a finite number of discrete 
charges give a better description. While the set of discrete charges are treated 
exactly, the continuum part is treated within a saddle point approximation to 
obtain a saddle point free energy, from which all the relevant transport quanti- 
ties are evaluated. The saddle point density of the continuum and its relation 
to the set of discrete charges show how the nature of the solution changes 
qualitatively with increasing disorder, and allows us to look for rare but al- 
lowed configurations that give rise to novel characteristics of the distribution. 
While two such characteristics, namely the existence of non-analyticity near 
g = 1 and the resulting asymmetry in the distribution have been reported 
earlier, we show here for the first time within the same framework that there 
exists a log-normal tail in P(g) for g <C 1 even for weak disorder where the 
average (g) >> 1, consistent with earlier field theory calculations. In order to 
check the accuracy of our method, we also obtain several other results in the 
metallic limit for which exact results are available. These include the average 
and the variance of conductance to order 1/r for all symmetry classes and the 
difference in the saddle point density between different symmetry classes. 

The method developed here is very general and can be applied to other quan- 
tities as well, e.g. the study of the distribution of shot noise as well as the con- 
ductance distribution in the presence of Andreev scattering in an NS (Normal 
metal-Superconductor) junction. 

One should keep in mind that the DMPK approach does not contain the effects 
of wave function correlations in the transverse direction, which are expected 
to be important in higher dimensions [28]. Nonetheless, the novel features in 
P(g) near g = 1 within our framework appears under very robust conditions 
which are independent of dimensionality. It is therefore possible that some of 
the features will persist in higher dimensions as well. The similarities of the 
shape of P(g) in the crossover regime obtained within the current framework 
with the numerically determined P(g) in 3D at the critical point appears to 
be consistent with the above conjecture. 
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A Partial densities 



In the following we give some details of the calculations for the partial densities 
mentioned in Sec. 4. 

It is convenient to define the integral: 

roo A 

n a (q)= dssm(qs)-V a (((s) + \ 2 ). (A.l) 
Jo as 

Then the partial density a a is given by 

°a = —\ dqcos(qt)(l-e-* q )n a (q). (A.2) 

• Partial density a z : 

Using the above definition with V z given by Eq. (46) we have 

^z(q) = / dssin(qs)- ~ (A.3) 

W Jo W ; dscosh 2 s + A 2 V ; 

which can be done by partial integration with result [29]: 

2 A(l + A 2 ) sinh (W2) 1 ' 



where 

a = arcosh(l + 2A 2 ) = 2arsinhy / A 2 = 2x 2 . (A. 5) 

Therefore the density a z can be written as 

a . (f) = i * n r dq s ^i sin c^i) e -» (a.6) 

w \/A 2 (l + A 2 ) 0s #7 1 V 2 > 

evaluated at 7 = n/2. From [29] and after algebraic simplification we obtain 
<j z as given in (59a). 
• Partial density a u \ 

We will first consider the metallic regime. Taking the derivative of V ul we 
have 

n ui = - 1 - r ds^^± (A.7) 
2 Jo sinh 2 s + A 2 -Ai v ; 
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and using sinh 2 s = (cosh(2s) — l)/2 then Q Ul is given by 

_ 7r cosh(ga /2) 

2 sinh ((7717 2) 

where A^ = A2 — Ai and a = n — 2arcsin y/Xj. Finally, after writing (1 — 
e~^ q ) = 2e -7n? / 2 sinh(7rg/2) in (A2), a u \ can be written as given in (59b) for 
A rf < 1. 

In order to obtain a u \ for A^ > 1, we follow a different method. First we 
take the derivative w.r.t. A^: 

dtt Ul 1 r°° , . , , d 



dssin(qs)— , 2 — . (A. 9) 



5Arf 2 Jo ds sinh s + A^ 

A partial integration gives the boundary term zero and 

dtt ui q dI Ql 



dX d 2 dX d 



(A.10) 



where 



dX d Jo cosh(2s) + {2X d - 1)' 



9IUl ~ ° r <** r7 TOS(9 ;\ v (A.H) 



(A.12) 



For X d > 1 [29], 

ggi = 7r sin[(g/2)arcosh(2A d - 1)] 
d\ d sinh(g7r/2) y / (2A d - l) 2 - 1 

Using arcosh(2Arf — 1) = y, we can integrate to obtain 

7T COS! QV /2 ) 

Jo, = 7-, ; — r + terms independent of X d . (A. 13) 

1 (7 sinh ((7717 2) 

Neglecting the A^-independent terms, then we have a result valid for A^ > 1. 
Therefore we have 



a. 



"Ml 



2sinh(g7r/2) 



cosh(ga /2), A d < 1; 
cos(g6 /2), A d > 1; 



where we have included the result for A^ < 1 obtained within this method 
(see Eq. (A8)) and b = arcosh(2Ad — 1). We can check that there is no A^ 
independent terms by comparing the two expressions of fl ui at the boundary 
X d = 1 i.e. at a = and b = for A^ < 1 and A^ > 1, respectively. We 
note that &o — ^ a 0i an d the expression for A^ < 1 becomes valid for A^ > 1 
if we allow a to be complex. The corresponding a Ul for all A^ is then given 
by (59b). 

Partial density a u2 
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Here we will use our approximation to V^ 2 (s), Eq. (58). This approxima- 
tion starts to fail for x 2 > 5 and/or ^ C 1, however the metallic limit 
x 2 — > is well described and also it has been shown in [2] that this approx- 
imation remains valid in the crossover regime. Equation (58) can also be 
written as 



KM 



i 



+ 



i 



2 

cts m s 



{s 2 + s 2 J 2 + Ac 2 s 2 



{s-ic) 2 + (3 2 {s + ic) 2 + f3 2 
with (3f = c 2 + s^. Introducing this expression in (Al), we get 
d r°° ds 



(A.15) 



0^(9) «-^/ o °° ^008(^)0) = -^ 

(A.16) 

where A± = (3± ± c ± as 2 n /2c, and the corresponding density becomes as 
given in (59c). 
Partial density ar 

In order to obtain ar we will use our approximation to V r (s), Eq. (57). 
In this case, we have 



fir~T / dssm(qs)s 
Jo 



1 - 



04 



ds 



cosh(2s) + cosh<5i 
cos(qs) 



d f 

= -n5 (q) + a 4 — / 

ag Jo cosh(2s) + coshoi 



—n5'(q) + a 4 ; 



7r (9 sin(g<5i/2) 



2 sinh <5i <9g sinh(g7r/2) 
Therefore ar is given by 



(A.17) 



2 z* 00 

<Tr(t) ~ — ^ / dqcosqt(l — e 
n 2 Jo 



-n5'(q) + 



7ra 4 9 sin(g5i/2) 



2 sinh 5 1 dq sinh(g7r/2) 

(A.18) 

The integrals for each term in (A22) can be done. The final result is given 
in (59d). 

Partial density a 2 

We first need to evaluate the potential V 2 given by (33), which can now be 
done with the densities evaluated above. Analytic expressions are possible 
only in the metallic regime, where the largest contributions come from ar 
(which is proportional to T > 1) and a z . Since V 2 itself is proportional to 
A2, we will ignore the contributions from a ui and a U2 as well as from 02, 
which will otherwise require a self consistent calculation. 

The potential V 2 r obtained from ar is 



v 2T = -rx 2 A(t) 



(A.19) 
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where we have defined 



A(t) = / 
J o 



ds 



t 2 



sinh(2t) sinh(2s) 



(A.20) 



and we have neglected the A2 dependent part of a-p- We extend the integral 
to —00 and continue to the complex upper half plane. Note that there is no 
pole on the real line. Using variable u = 2s and keeping in mind that poles 
on the real line are to be excluded, we can rewrite it as 



2Jc 



At 2 — u 2 sinh(w) 



(A.21) 



The poles are at u — ikir for k =integers. This gives 

k cos(kn) 



A = V 

^ k 2 + At 2 /n 2 ' 



(A.22) 



The corresponding f2 2 r, after a partial integration, can be written as 



7T 



n 2r = TX 2 — g [tanh(g7r/4) - 1] . 

o 

The density for this part is, after algebraic simplification, 



a 2 r = 



t 2 -n 2 



t 2 - 7T 2 /4 
it 2 + 7T 2 /4) 2 {f + 7T 2 ) 2 



(A.23) 



(A.24) 



B Partial free energies 

With the above partial densities of Appendix A, it is possible in principle 
calculate all the partial free energies F a b through the relation (48). However, 
only F zz , F ZU1 and F U1U1 are evaluated analytically exactly. For some of the 
remaining free energy terms we made some approximations which gives both 
the exact metallic and insulating regimes and a reasonable approximation 
(as checked numerically) in the intermediate regime, while other free en- 
ergies require different sets of approximations in different regimes. In the 
following we will point out the approximations used to obtain each expres- 
sion in the paper. 
Free energy F zz 

This can be obtained exactly. Although we have evaluated a z already, it is 
useful to start with the full expression of F zz in terms of the three integrals 
(Eqs. 48, 49 ) because as we will see, changing the order of integrations 
(e.g. doing the q integral last) will sometimes be more convenient. Using 
the result for Q z in (A4), writing cosh 2 t = (cosh(2t) + l)/2, and 1 - e _7r9 = 
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2 e -*5/2 s inh(7rg/2), we have 



F = 

- 1 zz 



2z 2 



Jx 2 (l + X 2 



dqqsm(qx 2 )e- nq/2 / dt 
Jo 



cos(gt) 



cosh(2t) 



+ 2A 2 
(B.l) 

The t integral can be done. Changing variables x = nq/2 and writing 
sin 2 (#) = (1 - cos(20))/2, the result is 



2tt(3 A 2 (l + A 2 ) Jo 







(2a \\ 


/ dx — 


1 — COS 


— X 


10 7T 




V 7T ). 



sinh a; 



(B.2) 



where a = 2x 2 has been defined for later convenience. The first term is 
known. The second integral can also be done by noting that (2x/ir) cos(2a7ra;) 
can be rewritten as {8/ da) sin(2a7ra;). The final result is given by (66), valid 
for all X 2 . This is a crucial term that dictates the crossover from metal- 
lic to insulating behavior. Note that the leading term for A 2 « 1, after 
cancellations, is 



F = 

A zz 



15(3 



[1 + 0(A 2 )]. 



(B.3) 



For A 2 = 0, this gives the exact variance in the zeroth order. 
• Free energy F^ z 

Using the expression for a z (Eq. 59a), we obtain by partial integration: 



Ft. = 



Yzx 2 



2a/A 2 (1 + A 2 ) 



f 

Jo 



dt 



ar sinh \J sinh 2 1 + A 2 

(cosh 2 1 + A 2 ) (sinh 2 1 + A 2 ) 
tsinh(2t) 
(t 2 + (7T 2 + a 2 )/4) 2 -a 2 t 2 ' 



x 



(B.4) 



We first obtain an expression valid in the metallic limit. 
Expand the integrand in powers of A 2 up to the first power. The leading 
A 2 -independent term leads to a simple integral, giving 



-,l(met) 
1 Yz 



zTxo 



zTx-2 



2y/X 2 (l + A 2 ) sinh(2x 2 ) ' 



(B.5) 



The next term involves 



dt 



( 1 



(t 2 + vr 2 /4) 2 sinh 2t 2{t 2 + tt 2 /4) 2 Vsinh 2 1 



(B.6) 



where a term 2(t 2 — 7r 2 /4)/(t 2 + 7r 2 /4) 4 in the integrand coming from differ- 
entiating the a = 2x 2 dependent factor in (B4) has been omitted because 
its integral over t is zero. The first term in (B6) is of the form 
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h=— 



d r°° dt 



dbJo t 2 + 6sinh(2t) 



d_ 

"db 



* (~l) fc 



7T^ 



1 + 2E 



'-V 



7T^ 



:[C(2)-1)]- 



(B.7) 
(B.8) 



Writing ^ 2 



d 1 



, ,,, , doing a partial integration and then 

using (f2+ f 2 /4)3 = (^2/4)2 [1 - (t 2+ 7r / 2 4 /4) ], the second term in (B6) can be 
rewritten as 

h = ~ 



dt 



sinh(2t) 



2t 



+ 



(t 2 + 7T 2 /4) 2 (t 2 + 7T 2 /4) 



,[3C(3)-C(2)-1], 



7T 



(B.9) 
(B.10) 



where C(x) is the Zeta function. Collecting terms, the free energy is 



F ( r<) w _^Txf _ f£x| [3C(3) _ 2C(2)]) (Rn) 



7T^ 



where we have used A 2 — > limit and kept only the leading x 2 dependent 
terms. 

To obtain an expression valid in the crossover regime, we use the fact that 
a z can be integrated exactly. We define 



= J dta z (t) 



2 sinh 2x2 



1 



1 



(t + X 2 ) 2 + 7T 2 /4 (t - X 2 ) 2 + 7T 2 /4 



Then 



F Vz = dtV r a z = V r cu z \™ - / dtV+u g . 
Jo Jo 



(B.12) 
(B.13) 



The first term is zero, and we use the approximation (57) for V£ to obtain, 
after simple integrals, 



-i(cross) 



sinh 2x2 2x 2 



where 
B = 



00 cosh z + cosh Si 



(B.14) 



(B.15) 



(Z + 2X 2 ) 2 + 7T 2 (^-2X 2 ) 2 + 7T 2 _ 

This can be done by contour integration. The integrand has poles at 

z = ±m±2x 2 - z = ±i(2ife + l)7r±£i, k = 0,1,2,- •• (B.16) 

Closing the contour on the upper half plane, the sum over residues can be 
rearranged and identified with DiGamma functions \l/(x): 
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cosh 2x2 — cosh 5\ 4x 2 — 5j sinh 5i 

+ ix 2 _] - - ix 2+ ] - + irr 2+ ] + - ^ 2 -]) , 



2 sinh <5i 



(B.17) 



where x 2 ± — x 2 ± Si/2. 
Free energy 

This is done exactly. Using the expression for V z in (46) and fl ui in (A14), 
we get 



F = 

- 1 ZUl 



dt 



7r Jo cosh t + X 2 Jo 



dqcos(qt) cosh(ga /2)e^ /2 . (B.18) 



We do the t-integral first, writing cosh 2 t in terms of cosh(2t). For the q- 
integral, we use e~ gn ^ 2 / sinh(g7r/2) = 2/(e qn — 1), giving 



F 



ZUl 



^2(1 + A 2 ) 



sinh a 



a 2 + a 2 , 2 cosh a — cos a 



(B.19) 



Finally, using Aj = sinh 2 Xj, a = 2x 2 , and cosh a — cos a = 2(1 + Ai), we get 



ZUl 



sinh(2x 2 ) 



4xo 



sinh(2:r 2 ) 



(B.20) 



4a; 2 + al 2 cosh x 1 _ 

As mentioned after (A 14), a is imaginary for Ad > 1. Using Q U1 for Ad > 1 
from (A14) gives the same result. 
The metallic limit of (B20) is simply 



pimet) 



ZUl 



z8x n 



:b.2v, 



Free energy F ZU2 : 

Using the expression for V z and fi U2 , writing cosh 2 t in terms of cosh(2t) 
and doing the t-integral first as in F ZUl results in elementary integrals over 
q. The result is 



TP ~ 

r zu 2 ~ 



2x 2 



4/?! sinh(2x 2 ) 



x\ + (tt/2 + /3i + c) : 



+ c -> -c 



(B.22) 



Note that f2„ 2 itself uses an approximate expression for V U2 . The free energy 
expression is a good approximation in both the metallic and the crossover 
regimes as checked numerically. In the metallic limit, A± — > j3\ — > Xd, c — > 
and (B22) gives 

^ (B.23) 



p(met) 

ZU2 



Free energy F z2 : 
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Since a 2 is valid only in the metallic limit, F z2 can be evaluated only in 
the metallic limit. We use the A 2 — > limit of V z to get 



-i(roet) 



zV\ 2 



z2 



^ C0s(/C7r) 



fc=l 



cosh t Jo 



dqqcos(qt) sinh(7rg/2)e I 2 9 ( 1+fc ). 

(B.24) 



Doing the ^-integral first and then the q integral, we get 

zTX 2 4 » cos(A;7r) 



(B.25) 



The sum can be rewritten as 

00 cos(/br) ~ (-l) fe 



E 



E 

A;=2 



-E- 

L k=l 



_ 1)k+ l 



k 3 



+ 1 



and evaluated in terms of Zeta functions, giving finally, for A 2 <C 1: 

^rA 2 



(B.26) 



p (wet) 
-^2 



7T^ 



■[3C(3)-4], 



(B.27) 



• Free energy Frr: 

This requires one approximation for the metallic limit, a second one for 
the insulating limit, and a third one for the crossover regime. 

For the metallic limit, we use expansions of Vr and of in powers of A 2 : 



-i(met) 

rr 



eft 



t 2 + A s 



2i 



sinh(2t) 



1 _ hi— ( I 



2 dt\t sinh(2t) t 2 + ^ / 

(B.28)' 

The first term diverges. However, it is independent of A 2 and so we ignore 
it. Rest of the integrals, keeping only terms up to order A 2 , gives 



rr 



16 



(B.29) 



where we have replaced A 2 <C 1 by x\. 

In the insulating limit, we approximate Vp* as 



Vr(«) 



r 

2 



^M^, for s « x 2 ; 

smh 2:2 
sinh(2s)s r „ 

— ^f, fors>:r 2 , 



Then we obtain or(£) for t x 2 and also for t ^> x 2 by dividing up 
the integral over s from zero to x 2 and then from x 2 to infinity, using the 
approximation for Vp\ We then obtain Fpp by dividing the integral over t 
into one from zero to x 2 and the other from x 2 to infinity. The final result 
is 



F, 



(ins) 



2^2 



rr 



pT'x 



(B.30) 
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In the crossover regime, we use (Jr(t) from (59d), and Vr(t) = Vp(t) + 
5V r (t), where K r ° = V T (x 2 -> 0) = Tt 2 /2. The first term in integral (48) 
coming from the first term of Up in (59e) involves 



Kl v = / dtV r (t) = / dt-Tt 2 + / dtSV T (t). (B.31) 
jo jo 2 Jo 

The first term diverges; but it is independent of Xi and so we ignore it. We 
do a partial integration for the second term and use the approximation (57) 
for 8V^, leading to 



Kl r = a 4 r / 

JO 



dt 



t 2 



cosh 2t + cosh 5i 



— CL4F 



1 5 r» 
A <9A 



r 

JO 



dt 



cos At 



cosh 2t + cosh Si 



J A=0 



(B.32) 

Expanding the result for the integral in power series in A and taking the 
limit we obtain 



a^T 81 , 2 , X 2\ 



A'. - — 

rr 24 sinned 

The second term of <xr leads to integrals of the form 

1 1 



it 



= / dtVAt) 

JO 



t 2 + c 2 t 2 + c 2 



(B.33) 



(B.34) 



The dominant term Vp(t) gives r7r(c_ — c + )/4. The next term 5Vr is again 
done by partial integration, leading to 



Ya A / dt- 

JO 



cosh 2t + cosh <5i 
The third term of <7r leads to 



1 t 1 t 
— arctan arctan — 

c + c+ c_ c_ 



K 3 
i\ rr 



1 



cosh 2£ + cosh <5i 
Again the dominant term of Vr(t) gives 

The next term, after a partial integral, leads to 
Ta 4 f°° t 



(B.35) 



(B.36) 



(B.37) 



sinh 8 



-f 

1 jo 



dt 



cosh 2t + cosh 81 



tanh 1 tanhttanh 



81 



(B.38) 



Collecting all the terms with their appropriate pre-factors, we get the ex- 
pression 



Frr — 



T 2 a A 8 l 
12 sinh 81 



(- 2 + si) + 



r 2 ai 



dt 



t*{t) 



in sinh 81 Jo cosh 2t + cosh 8 1 



, (B.39) 
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where 



_ . . 7r + iSi 2t 
9(t) = — - — arctan 



71 — 15\ 2t 

■ arctan ■ 



n — i5\ 2 7r + iSi 

— 7r tanfT 1 [tanh t tanh(5x /2)] 



(B.40) 



The integral with can not be done. However, to a good approximation, 
we can write 



*(*) 



£>if 



(1 + sflhlDd 2 ) 



Am5i . Si .tt^i 2 2 

" ; Dl= 2 i x2 ~ mtann Tr; ^2 = (vr 
! ir z + di 2 12 

(B.41) 



Denoting 77' = yDi/D 2 , we need the integral 
1 Dxt 2 



/ 



dt- 



D 1 d 
'^2~dr) 



cosh 2t + cosh 5x (1 +r/t 2 ) 2 
This integral is similar to (B15), and the result is 
dt 1 1 



dt 



cosh 2t + cosh 5i 1 + r/t 2 
(B.42) 



00 cosh 2i + cosh Si t 2 + a 2 ia sinh 5i 

1 a i5\ 



x 



I 2 71 271 



1 + a_iS 1 \l 

2 71 2n h 



Note that the integral is true for Si either real or imaginary. Collecting 
terms, we finally have 



Frr — 



(- 2 + si) 



X 



where 



T 2 a 4 Si 
12 sinh Si 



1 77 ± i<Ji/2 



0.44) 



(B.45) 



and is the derivative of the DiGamma function. 
Free energy F Fui : 

This energy is analytically calculated in the metallic and insulating regimes 
and numerically in the crossover region. 

In the metallic limit, we use a Ul and the A 2 — > limit of Vp: 



-i(met) 



Tu-i 



An Jo [t 2 + fl 2 _ t 2 + fi 2 + _ 



7\ ± ao 



(B.46) 



Rewriting /ijt 2 / (t 2 + jif) — ^ (1 — /i 2 /(t 2 + /i 2 )) and the relation //_ + /i + = 



43 



7i, we get 



{me t)^JT [<*> 



( /£ n% 

- T7T- <5~ + 



(B.47) 



The first term is divergent but constant, so we ignore it. The second term, 
after using the definitions of //j, gives 

Ft £t) - f 4 (B-48) 

Using a = 7r — 2 arcsin ^/A^ and <C 1, and ignoring constants we finally 
obtain 

- "X^- (B.49) 
In the insulating limit, a"f m ^ w 2r/:r 2 which gives 

^1T } - -^3x 2 . (B.50) 

Free energy F ru , 2 : 

Fy U2 is calculated only in the metallic limit while for the crossover region, 
it is obtained numerically. Contribution in the insulating limit is negligible 
compared to F^™\ and will not be evaluated. 

In the metallic limit we proceed as for F rui . Using a U2 , the dominant 
term V®(t) = Tt 2 /2 (defined above (B31)) gives a diverging constant that 
we ignore, plus a term 

.(met) ^ Pr^Pl 

4 

where we have kept only the leading term in powers of A 2 <C 1. Since in this 
limit Pi Xd, we obtain 



F (met) ^ (B 51) 



^r« 2 ~ -^TT^d. (B.52) 



Free energy F U1U1 : 

This can be done exactly. We use V Ul from (46) and Q Ul from (A 14). We 
do the t integral first, 

PCX) 

Aiiui — / c?tln(sinh 2 t + Ad) cos (<?£). (B.53) 

JO 

A partial integration gives a diverging term. We introduce a cut off b: 

t 2 ih ■ i uw 7r cosh(ga /2) 

iuiui = - psm gfi k-^oo — — — B.54 

q q smh(g7r/2) 

The g-integral for the free energy leads to a diverging constant from the 
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first term and we ignore it. The second term can be rewritten in the form 



P f°°dq _ m?/2 l + sinh 2 (ga /2) 



p f 00 dq 
~2 Jo ~~q 



q sinh(g7r/2) 
Using x = qn/2, this integral can be done, and we get 



(B.55) 



mm 



I 
2 



iln 

2 



a 



sinao 



+ 



f 

Jo 



dx e~ 



x sinhx 



Neglecting the diverging constant, we finally have 



a 



sin a 



(B.56) 



(B.57) 



Again, a is imaginary for A^ > 1. 
Free energy F U2U2 : 

Given the form (59c) for a U2 , this can be done exactly. However, a U2 
itself has approximations. The resultant approximate F U2U2 will be valid in 
the crossover regime as checked numerically, and will also give the correct 
metallic limit. 

Using <y U2 and a partial integration, we start with 



P 



i lo 



dtVL x 



A ( t t 

A + arctan arctan — 

V Ai A 2 



A ( 1 1 

+ A_ arctan arctan — 

V A 3 A 4 



,(B.58) 



where A 1 — fix + c, A 2 = ir + P± + c, A 3 = fix — c, A 4 = 7r + Pi — c. Consider 
the term with A^ Derivative w.r.t. Ai gives 



9Ai U2U2 



MPi 



dtVL 



t 



l J—oo 



t 2 + A\ 



(B.59) 



The integral can be done by simple contour integration. Then integrating 
back w.r.t. Ai gives the final result 



U2U2 SttPi 



where we defined 

7T 



In 



(A i + p 1 y- c 2 



+ 



7T 



C — 



CL\S n 



In 



K + Pi + c 
Ai + P ± -c 



(B.60) 



(B.61) 



In terms of these definitions, collecting similar terms from other A, terms, 
we finally have 



«2«2 



P 

Mp^ 



[A + (h Al - h A2 ) + A_(h As - h A J] . 



(B.62) 
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In the limit A2 <C 1, A± — > — > 7rln[Aj + and /?i — ► Xd, so that 
the metallic limit of F U2U2 is given by 

F<££>**£]nx d , (B.63) 

where we have omitted the constant terms. 
Free energy F UlU2 : 

This again uses the approximate a U2 , which is valid in the crossover regime 
as well as the metallic limit. We proceed like F UlUl . Using V Ul and Q U2 , and 
doing the t- integral first, we generate a diverging constant term that we 
ignore. The rest becomes 

F U1U2 » - A r ^ cosh(ga /2)e-^/ 2 )[A + e^ A + + A_ e -« A -] (B.64) 
4pi jo q 

where A± — fli ± a Taking a derivative w.r.t. ao we get rid of the 1/g term 
and replace the cosh(gao/2) factor with sinh(gao/2). The integral can then 
be done, which can be integrated back w.r.t. ao, giving 



F ~ A 

U1U2 ~ 8A 



A + ln[(7r + 2A + ) 2 -a 2 ] + A_ln[(7r + 2A_) 2 -a 2 ] . (B.65) 



As before, for > 1, ao is imaginary. In the limit A2 <C 1, we get 

FiTut^^d. (B.66) 

Free energy F 2 r' 

Since V 2 is already only up to order A2, this term is only valid in the 
metallic limit. We use the A2 — > limit for ar- 

TB 

F 2r = -f (V 2r + V 2z ) . (B.67) 

The first term F^ r involves an integral over A(t) defined in (A20) which can 
be easily done, and using Y^=i cos(kn) = —1/2, gives 

F r (met) w r 2 ^ (B 6g) 



The second term involves several parts. The first one is 

roc q 2 _ 3^ 2 

■ h= L WT^Y> Am ' = I/2 ' 



7T 2 

= k cos(A;7r) — 
k=i 2 



evaluated at 6 = 7r 2 /4 and c = n 2 k 2 /A. The result for the integral in (B69) is 
(4/7T 2 ) YjT=i cos(k7r)/(k + l) 3 , after taking the derivatives and using values 
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of b and c. The sum can be obtained in terms of Zeta function, giving 



J 1 = ^[3C(3)-4]. 



7T 



The second part is the same as in (B8), given by 

J 2 = ^[C(2)-1)]. 



The third part is of the form 
dt t 



f 

Jo 



l d 2 dt t 



(t 2 + a 2 ) 3 sinh(2t) 2 db 2 Jo t 2 + b sinh(2t) ' 



/ 

Jo 



(B.70) 



(B.71) 



(B.72) 



which can be done using the previous results. After the derivatives it be- 
comes 



'-I 



-1) 



3 + 4 £(ETi? + 2 £(ETi? 



-[3C(3)+C(2)-3)]. 

7T^ 



The other integrals are elementary, of the form 



Jo 



00 dt 



(-1) 



r-1 



7T a 



l-2r 



7T 



{t 2 + a 2 Y v ^ 2 (r- l)!r(3/2-r)' 
Collecting all the terms, we finally have 

Thus the Free energy contribution is 



(B.73) 
(B.74) 



^B.75) 



-iz(met) 



[4-3C(3)] 



(B.76) 



C correction to the functional integral 



Here we give some details for the deduction of the relation (65). 

As we mentioned, in general the integral (64) is of the form I a ~ 1/ \Jdet\u(t, s)\. 
In the regime A2 <C 1, we can expand u as u = uq + A2<5w, where uo is inde- 
pendent of A2. Therefore the determinant can be also expanded as det\uo\ x 
det[l + A2M0 The first term is independent of A2, and only contributes to 
the normalization. The second part can be written as exp[Tr ln(l + XzU^Su)] 
Expanding the logarithm, it is clear that the correction to the Free energy is 
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going to be proportional to A2, and the proportionality constant is small be- 
cause the dominant metallic contributions come from the Uq part and 5u is a 
small correction. Since all other terms in the Free energy in the metallic limit 
are of order T-^X^ ~ 1, the fluctuation corrections for the functional integral 
are in general negligible. 

However, our saddle point density goes to zero at the lower limit t — 0. It is 
approximately constant, equal to T, for large t. The scale over which it changes 
is t 2 ~ (A 2 — Ai). In general, for the saddle point value of a less than the width 
of the peak in u vs a, the lower limit of a(t) = will matter. For example, 
in the extreme case where the saddle point a — 0, the integral over a would 
only give 1/2 the full value. We can not evaluate this correction explicitly, but 
we only need the Ai and A2 dependent contributions, which we can estimate. 
Let us consider the region of integral over t, s up to some constant a times 
\ l J 2 = (A 2 — Ax) 1 / 2 , which is the scale beyond which the t dependence becomes 
negligible. This is the part where the correction comes from. The integral I a 
will have a contribution coming from the exponent 

.1/2 .1/2 

J a = dt dsSa(t)Sa(s)u(t,s). (C.l) 

JO JO 

1 /2 

We scale out the A^ from the limit by changing variables t' = t/a\ d , s' = 

1/2 

s/a\ d , which leads to J a = \d,J' a where we defined 

J' = a 2 f 1 dt' f 1 ds'5a(t')5a(s')u(t', s'). (C.2) 
J J 

Since u is logarithmic in t, s for small t, s, the change of variables will produce 
a In Arf term, plus terms independent of A. Expand 

00 00 
5a(t)= £ e^ nt a n , u(t,s)= £ e^ mt+ns \ mn , (C.3) 

n=— 00 m,n=— 00 

such that J' a ~ J2-oo VmUmnCTn, where 

u mn = C ds' f 1 dt'e i2nims ' +nt,) u(s', t'). (C.4) 
JO Jo 

It is now clear that if u oc In \f\~d~ is the dominant part, then only m = 0, n = 
part is non-zero because the integral is proportional to terms like sin(27rm)/m. 
In other words, u mn = uoS m fiS n fi, where Uo = — In y/Xj. Then we can write 

since a sp 0. (The correction will again be powers of A). In the Free energy, 
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this gives a contribution (Eq. 65) 

5F « -lnA d . 
2 

Note that because of the logarithm, the coefficient of the In term is one-half, 
despite all the approximations. The term appears simply from the scaling of 
the original variable, which was determined by the form of the diverging part 
of o~ UljU2 which sets the scale over which the saddle point density becomes 
independent of t near the origin. Note also that for ^> 0, the density no 
longer has a dip at t — 0, and the correction vanishes. Thus it will be relevant 
only in the metallic limit. 
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